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Abstract 

In this paper we consider the question of the degree to which negative and positive 
energy are intertwined. We examine in more detail a previously studied quantum state of 
the massless minimally coupled scalar field, which we call a "Heifer state" . This is a state in 
which the energy density can be made arbitrarily negative over an arbitrarily large region 
of space, but only at one instant in time. In the Heifer state, the negative energy density 
is accompanied by rapidly time-varying energy fluxes. It is the latter feature which allows 
the quantum inequalities, bounds which restrict the magnitude and duration of negative 
energy, to hold for this class of states. An observer who initially passes through the negative 
energy region will quickly encounter fluxes of positive energy which subsequently enter the 
region. We examine in detail the correlation between the energy density and flux in the 
Heifer state in terms of their expectation values. We then study the correlation function 
between energy density and flux in the Minkowski vacuum state, for a massless minimally 
coupled scalar field in both two and four dimensions. In this latter analysis we examine 
correlation functions rather than expectation values. Remarkably, we see qualitatively 
similar behavior to that in the Heifer state. More specifically, an initial negative energy 
vacuum fluctuation in some region of space is correlated with a subsequent flux fluctuation 
of positive energy into the region. We speculate that the mechanism which ensures that the 
quantum inequalities hold in the Heifer state, as well as in other quantum states associated 
with negative energy, is, at least in some sense, already "encoded" in the fluctuations of 
the vacuum. 

PACS numbers: 04.62. +v, 42.50.Dv, 03.70.+k, ll.10.-z 
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I. INTRODUCTION 



It has been known for quite some time that quantum field theories generically 
contain states associated with negative energy [if. Indeed, negative energy seems to 
be required by the laws of physics for such effects as Hawking evaporation (which 
allows the unification of black holes with the laws of thermodynamics) 0] , and for the 
stability of the Minkowski vacuum. In the former case, the positive Hawking radiation 
observed at infinity is paid for by a flux of negative energy down the horizon of the 
black hole In the latter case, since the vacuum is not an eigenstate of energy 
density, there will be energy density fluctuations in the vacuum. In order for the 
averaged energy density to be zero, there must exist negative, as well as positive, 
energy density vacuum fluctuations. States associated with negative energy are now 
routinejyproduced in the laboratory, e.g, squeezed vacuum states and the Casimir 
effect |4|, |5[ . However, the energy density itself is too small to be directly measurable. 

On the other hand, unrestricted amounts of negative energy could result in gross 
macroscopic effects such as violation of the second law of thermodynamics |6 , [7 1 , and 
exotic spacetime geometries such as wormholes (1, 0], warp drives [l(J II, U 13], and 



time machines [14J. However, the laws of quantum field theory contain restrictions 
on negative energy in the form of "quantum inequalities" . These involve bounds on 
the magnitude of negative energy and on its temporal or spatial extent. There has 
been much progress in this area in the last 15 years (for some recent reviews see 



Refs. 15, 16 



Some years ago, Heifer [18| suggested the existence of a class of quantum states 
in which one could make the energy density arbitrarily negative in an arbitrarily 
large region of space, even in Minkowski spacetime. We call these "Heifer states." 
The present authors, together with Heifer, verified this claim in an earlier paper 
(l9| . which we will refer to as FHR. In that paper the characteristics of the energy 
density were analyzed. We argued that a crucial feature of these states is the fact that 
although the energy density can be made arbitrarily negative over an arbitrarily large 
region of space, this can be done only at one instant of time. The worldline quantum 
inequalities must hold for these states as well, since they have been proven to hold for 
all quantum states. Therefore, for this to be true, there must be rapidly time-varying 
fluxes which accompany the negative energy density. In the current paper, we show 
that this is indeed the case. The expectation values of the energy density and flux 
in a particular Heifer state are calculated and graphed, and their correlations are 
studied. An inertial observer who initially passes through the negative energy region 
must quickly encounter fluxes of positive energy which subsequently enter the region, 
thus ensuring that the time-averaged sampled energy density along the observer's 
worldline is non-negative, in accordance with the quantum inequalities. 

These fascinating correlations of flux and energy density are by no means unique 
to the Heifer states. In the second part of the paper we calculate the energy density- 
flux correlation function for the Minkowski vacuum state, in both two and four- 
dimensional spacetime. Although here we analyze correlation functions, as opposed 
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to expectation values, nonetheless, we find remarkably similar behavior for the fluc- 
tuations of energy density and flux in the vacuum to their expectation values in the 
Heifer states. An initial negative energy vacuum fluctuation in some region of space 
is correlated with a subsequent flux fluctuation of positive energy into the region. We 
speculate that the mechanism which ensures that the quantum inequalities hold in 
the Heifer state, as well as in other quantum states associated with negative energy, 
is, at least in some sense, already "encoded" in the fluctuations of the vacuum. 



II. THE HELFER STATE 



Here we briefly summarize the results of FHR, and refer the reader to that paper 
for more details. The quantum state used in FHR is a superposition of two-particle 
states, which can be expressed as 



N 



+ 



J d 3 h d 3 A; 2 6(k 1 ,k 2 )|k 1 ,k 2 ) 



where the two-particle state, |k 1? k 2 ), contains particles with three-momenta kx and 
k 2 . The normalization factor N is given by 



N 



We consider the case where 



1 + 2 / (PhSk^b^M) 



6(k 1 ,k 2 )= X (k 1 + k 2 )(|k 1 ||k 2 |) 



-i 



and 



Xo, if |p| < Po 
, otherwise 



(2) 



(3) 



(4) 



where xo and po are arbitrary constants, and p = ki + k 2 . (This is the v = —1/2 limit 
of the state discussed in FHR.) In addition, we require that the magnitudes of the 
momenta of the particles be less than a cutoff parameter, A, so that |ki|, |k 2 | < A. 
Thus, the state is described by the three parameters, A, xo, an d Po- In the limit that 
A becomes large, with xo an d Po fixed, the two members of the pair of particles have 
nearly opposite momenta. This is the limit of particular interest. 



A. Energy Density in the Heifer State 

The energy density of a massless, minimally coupled scalar field may be obtained 
from the renormalized two-point function using 

p= l - lim [{d t d t , + V • V) (: <p(x) <p(x') :)] . (5) 
t'—>t 
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Units in which % = c = 1 are used throughout this paper. In our case, this two-point 
function may be shown to be 

27V 2 r 1 
(: <p{x)<p{d) :) = 7^ Re / d3kd3k ' 77=^1 



LOUO' 



x 2e i(k ' x '- k - x) e i(w '-"' < ' ) /d 3 A;i6*(k 1 ,k)6(k 1 ,kO 

+ e i(k.x+k'.x') e -i(o, (W t') fe ( k) k /) j _ (g) 

The resulting energy density is 

p = ^pRe J ftft'V^/(l + k.k') |- 2e i(k'-k).x e i( w - w )t y d3fc 1 &*(k 1 ,k)6(ki,k') 
_ e «(k+k').x e -<( w w)t 6 ( kj k ')j ; ( 7 ) 

where k and k' are unit vectors in the directions of k and k', respectively. It is 
convenient to write p = pi + P2, where p\ is the contribution quadratic in b, and P2 
is that linear in b. 

In general, it is difficult to evaluate p explicitly for the choice of &(k 1; k 2 ) given by 
Eqs. (13]) and (jl]). However, in FHR approximate forms were derived in the limit that 
the dominant contribution comes from modes whose frequency is large compared to 
Po, which is expected to be the case when A ^> p . These approximate forms contain 
an undetermined constant q ^> 1 which is introduced by making a lower cutoff in the 
a;- integrals of qpo- The result for p 2 is 

P 2 ~ ~^§~2~ fziPoir) ! diOLd' 1 cos(2ujt) , (8) 

07T Jqp 

where 

Air 

fi(jPo, r ) = [ 3 (Po r2 - 2 ) sin(p r) - p r(plr 2 - 6) cos(p r)} • (9) 
The function f2{po,f) has the following behavior: 

\ - (An/ r 2 )p 3 cos p r, r > pcT 1 
h(Po,r) ~ < ■ (10) 

[(47r/5)pg, r^po- 1 

At t = 0, Eq. © becomes 

— — / 2 (po,r) In 

D7T 2 V^O/ 

Let us now examine the behavior of p\. From Eq. (33) of FHR, we have that 



P2«-^ =r /2(P0,r)ln( — J . (11) 



Pl = 2 *f N2 Re / d 3 kl \ kl \-* I S p I d 3 p > (i + k ■ k' 



(27r) 3 J Ap\<po Ap'\<po 

p _ kl |- 1/2 | P ' - kl \-^ e -i(p-p'>x e il(\p-^\-\p'-*i\)t] > ( i2) 
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where p = ki + k and p' = ki + k'. In the high frequency limit, p\ at t = is given 
by 

Pi-^^A(P0,r)ln(A) , (13 ) 
tH \qp J 



where 



and 



/i(Po, r ) = — — [sin(p r) - Por cos{p r)} 2 , (14) 



(4vrp ) /r 4 cos 2 (p r) , r > p , ir , 
/i(Po,0 ~ ^ . (15) 

(4vr/3) 2 po, r^po" 1 

The time dependence of p\ is given by the factor of e^lP-^Hp'-kil)*] j n Eq. flT2J. In 
the high frequency limit, 



|p - ki| = ki 



^1-2^ + ^^-p.k,, (16) 



where fci = |kjj. As a result, 

u - J = | p _ k x | - |p' - ki| « ki • (p' - p) . (17) 

Because |p|, |p'| < po, the time dependence of p\ is on a time scale of order 1/po, and 
can be neglected if we are interested only in times for which \t\ -C 1/po- The time 
variation of the regions of negative energy is governed by the rapid time dependence 
of p2 coming from the cos(2u;t) factor in Eq. (jSJ). This is the time variation which 
allows the quantum inequalities to be satisfied in this case. 
The integral in the normalization factor, Eq. (j2J), is 

I = xt[d 3 k[ d 3 p(|k||k-p|)- 2 , (18) 

where k = ki and p = kx + k 2 . In the limit where |k]J ^> Pq, (|k||k — p|)~ 2 ~ |k|~ 4 , 
so we can write 



I « Xo / d 3 k |k|- 4 / d 3 p = Airxl ^ u- 2 du [ d 3 p 
When A 3> qpo, this becomes 



16tt 2 2 /po 2 _ po^ 

3 Xo \q A 

(19) 



167r 2 9 p n 2 



which leads to 



Nx{l + ^Zf'\ (21) 
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p 




-0 .0005 



FIG. 1: The total energy density p = p± + P2, is plotted as a function of r, at t = 0, 
for varying A. We have chosen po = 1 and q = 10. The energy density becomes more 
negative with increasing A. Note that on scales where r < Po 1 = 1, the energy density is 
approximately constant and negative. Such a region may be made as large as one likes by 
decreasing pq. 

For any value of po, the energy density is approximately constant in space over a 
region of size, r < Po~ l - In this region, (/2//1) ~ 9/(20 vrp ). If we choose 

80 7rp 

then IP2I/P1 > 1 and hence p < in this region. In the limit that A — ► 00, we can 
make p arbitrarily negative in this region. Furthermore, by making po small, we can 
make this region arbitrarily large. Note that if Eq. (|22|) is satisfied, and if q ^> 1, 
then N pa 1. In this case, the unknown parameter q and the cutoff parameter A only 
appear in the energy density in the ratio A/q. Thus we can now take this ratio to be 
the effective cutoff. 

For the purpose of plotting the energy density, we let q = 10 and set 

3 . . 

Ao = ™ , 23 

800 npo 

which ensures that there is a region of negative energy about r = 0. In Fig. [TJ we 
show the energy density as a function of r at t = for several values of A. In Fig. [21 
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FIG. 2: The total energy density p = p\ + P2, as a function of r, for varying t. We have 
chosen A = 1000, po = 1 and q = 10. At t = the energy density in the region near r = 
is negative. By t = 0.005, it has become substantially less negative, and by t = 0.05 it has 
already become positive. 

we illustrate the energy density for one choice of A at several values of t as functions 
of r. The special feature of this quantum state is that it exhibits large negative 
energy density at t — 0, but this negative energy rapidly disappears, as required by 
the quantum inequalities. The time scale for this change is small compared to the 
light travel time across the negative energy region. This will be discussed in more 
detail in the next subsection, where we show that this rapid change in p is associated 
with large energy fluxes. 



B. The Energy Flux in the Heifer State 

Let us temporarily switch to Cartesian coordinates. Then in the z-direction we 
have 

(T 0z ) = — (T° 2 ) = 1 lim [(d t d z , + d t ,d z ) (: <p(x) V {x') :)] . (24) 
Using Eq. ([6]), we obtain 

(T 0z ) = h + I 2 , (25) 
8 



where 



N 2 



(2nf 



3 f {Uk' z + U'k z ) J{ 



LOUO' 



and 



I 2 = Re / d A k J d 6 k x = c " • v 

2N 2 



k'+k).x e -i(uM-a/)t ^ k /j ? ( 2g ) 



h 



Re / cf / d A k 



,zu (uk' z + Lv'k z ) nil 



k'-k)x i(w-ui')t 



(2tt) 3 

x y"rf 3 A; 1 6(k 1 ,k)6(k 1 ,k') • 



LOUO' 



(27) 



In the last expression, we use the fact that 6(k, k') is real. At t = 0, we have 



h(t = 0) 



2N 2 



d A k \ d A k 



3 (ujk' z +uj'k z ) 



(27r) 3 7 "J v^Cc" 

x / rf 3 fc 1 6(k 1 ,k)6(k 1 ,k'). 



cos[(k' — k) ■ x] 



(28) 



However, = 0) = 0, as may be seen by making the replacements k — > — k, 
k' — ► — k', and kx — > — k^ and using the fact that 6(k 1; k) = £>(— k 1; — k). We can see 
from Eq. (ITT)) that 7i will only vary on a time scale of order 1/po, as does p\. Thus, 
if we are interested in times close to t — 0, we may set I\ ~ 0. 

Next we proceed as in FHR and examine h- Using p = k + k', and taking the 
limit in which to becomes large compared to po, we can expand the prefactor in h as 
follows: 



(uk' z + Lo'k z 



LOLO 



k| |p-k| 



(p-k)-z 



p-k).(p-k) 



kz 



LO 



k| |p-k| 

» p* - (k • p)fc 2 + ■ ■ • , 
where A; z = k • z. So we can now write the integral h as 



LO 



LO' 



(2tt)' 



Re / d 6 k / <f p [p* - (k • p)jy e jp x cj" 2 e 



2 -2uot 



|p|<Po 



Now let 
where 



ha + -^26 • 



ha = 



XoN 2 
(2tt) 3 



Re / d 3 k f d 3 pp z e ip - x to- 2 e~ 2iuJt 



(29) 

(30) 

(31) 
(32) 
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and 

I 2b = Re / d 3 k [ d 3 p (k ■ p)* a e^ x ^ 2 e" 2 ^ . (33) 

(2-k) 6 J J\p\<po 

First let us evaluate J2 a . Performing the angular integrations, we have 

d 3 k = 4vr / duu 2 (34) 



and 

i*,EE / d 3 p Pz e ip - x = f d 3 p Pz e iprcose , (35) 

where r = |x| and # is the angle between x and p. 

We are assuming spherical symmetry, so let us choose, without loss of generality, 
x to point along the z-axis. Then p z = p • z = pcos8, where p = |p|. We can now 
write 

Aa= I d 3 pp z e iprcose =2n f° dpp 3 P dcce iprc = -i g 2 (p ,r) , (36) 

J\p\<Po JO J-l 

where c = cos 9 and 

02 (Po, r ) = ^1 [ 3 Por cos(po^) + (poV - 3) sm(p r)] . (37) 
Therefore /2a becomes 

ha = -j^xg 2 (po,r) [ du sin(2c^) , (38) 

(z7Tj 7gp 

where as before, g is a constant chosen so that the approximations made in finding 
the integrand (u 3> po) are valid throughout the range of integration. Evaluating the 
integral we have 



T X ° N2 ( \ 

ha = -^ 2 -92{Po,r) 



cos(2At) — cos(2gp ^) 
t 



(39) 



We now turn our attention to hb- First fix the direction of p and integrate over 
the directions of k. In a right-handed x, y, z coordinate system, define: 9' to be 
the angle between k and the z-axis, 9 to be the angle between p and the z-axis, <$' 
to be the azimuthal angle in the x, y plane for k, <ft to be the azimuthal angle in 
the x, y plane for p, and 7 to be the angle between k and p. Then we have that 
k z = k ■ z = cos#',k ■ p = |p| COS7, and p ■ x = prcos9. Making use of the fact 
that 

COS7 = cos# cos9' + sin9 sin#' cos(y? — tp') , (40) 
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we can write 



1-21, 



XoN 2 
(2tt) 3 

Xo^V 2 
(2tt) 3 



Re / d 3 p e ip ' x / d 3 kpcosj cos9' 

J\p\<Po J 



w -2 e -2iart 



Re 



d 3 p e ipx ^p cos# I due 



Air 



\p\<Po 



(27T) 



Re 



Therefore our result for .Z2 is 



p|<Po 



-2iu)t 



d 3 p p z e ip ' x due 



(41) 



-^2 — ha + ^26 — g^2a ; 



and so 



67T 2 



cos(2At) — cos(2gp ^) 



(42) 



(43) 



where ^(po; r ) is given by Eq. ( 1371) . 

It is instructive to plot the energy density and flux at constant r as a function of 
time. The flux F is given by 



F=(T ^> = -(T 0z > = -/ 2 . 



(44) 



In Figs. [3] and HJ we plot the energy density and flux at fixed r = 2 as a function of 
time. In both plots, we have taken po — 1 and q = 10. Figure [3] has A = 100, while 
Fig. H]uses A = 1000, resulting in a more rapid variation in time. 

In Fig. [31 as we proceed from negative to positive values of t, we see that the energy 
density starts out approximately constant and positive. It then dips and eventually 
becomes negative, and afterwards becomes positive and constant again at large t. 
Note that the dip in the energy density is accompanied by a large positive outward 
flux away from r = 2, or equivalently a large negative inward flux towards r = 2. The 
flux goes through zero when the energy density reaches its minimum value at t — 0. It 
then becomes a large and negative outward flux, or equivalently a large and positive 
inward flux towards r = 2 as the energy density rises and becomes positive again. 
Figure H] shows the more rapid oscillations of the flux when A is increased. Note 
that in this case, the energy density undergoes small oscillations superposed upon 
its approach to a large negative value at t — 0. As originally discussed in FHR, the 
energy density may be made arbitrarily negative over an arbitrarily large region at 
one instant of time. We see here that the accompanying rapidly oscillating fluxes will 
ensure that the quantum inequalities are not violated, by quickly "pumping" enough 
positive energy into the region such that the time integral of the energy density over 
any observer's worldline will always be positive. 
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FIG. 3: The energy density, p, and flux, F, at r = 2 are plotted as functions of time. Here 
we have chosen pq = 1, q = 10, and A = 100. In order to be visible on this graph, the energy 
density has been multiplied by a factor of 50. Our plot is over a timescale small compared 
to l/po = 1, hence we can ignore the long-timescale dependence of F and p. Note that the 
sign of the flux changes when the energy density reaches its maximum negative value. 

Recall that the spatial region over which the negative energy density appears can 
be arbitrarily large. In particular, the light travel time over this region can be large 
compared to the time scale on which the energy density changes sign. One can see 
this in our plots by comparing Fig. [1] with Fig. HI The time scale for the variation 
of the energy density is about two orders of magnitude smaller that the size of the 
negative energy region. Thus one cannot think of the flow of positive energy into the 
region as beginning at the boundary and moving inward in a causal way. Rather, 
the energy density and flux are correlated over spacelike intervals. The energy fluxes 
depicted in Figs. [3] and H] are coherent over the entire region. A related example of 
spacelike correlations in quantum states with negative energy density was discussed 
in Ref. 0. 
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FIG. 4: The energy density and flux at r = 2 are plotted as functions of time, with po = 1, 
q = 10, and A = 1000. In this graph the energy density is multiplied by a factor of 300 in 
order to be easily visible. Our plot is over the same timescale as in Fig. [3j The increase 
in A results in a higher frequency of oscillation for the flux, F. This in turn results in the 
very short timescale fluctuations of the energy density seen in the figure. 

III. ENERGY-FLUX CORRELATIONS IN THE VACUUM 

In this section we discuss a different but related issue - the energy density-flux 
correlations in the Minkowski vacuum state. We saw in the previous section that, 
in the Heifer state, the regions of negative energy on a t =constant surface were 
correlated with a rapidly time- varying energy flux. Indeed it is this correlation which 
ensures that the worldline quantum inequalities are satisfied. As the energy density in 
a spatial region became negative, there was an accompanying positive energy outflow 
from (or equivalently, a negative energy inflow into) the region. When the energy 
density in the region began to rise and become positive again, this was correlated 
with a positive energy inflow to (or equivalently, a negative energy outflow from) the 
region. 

Does this behavior represent the general case? One might think that the answer is 
no. The emission of spatially and temporally separated pulses of ne gati ve and positive 



energy by moving mirrors would not appear to fit this profile [22|, |23|, |24(. However 



the moving mirror example is two-dimensional. It has been argued elsewhere that, in 
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four-dimensions, the emission of two separated plane sheets of negative and positive 
energy would violate the quantum inequalities, unless the positive sheet were infinite 



in extent [21]. So the spatial isolation of negative and positive energy in the moving 
mirror case may well be an artifact of two- dimensions [251 ] . 

So far all of the scenarios we have discussed involved expectation values of the 
energy density and flux. In this section, we analyze energy density-flux correlation 
functions in the Minkowski vacuum state, in both two and four dimensions. We find 
an analogous correlation between energy density and flux here as well. That is, if an 
energy density fluctuation in a region of the vacuum is initially negative, it is likely 
to be followed by a flux fluctuation of positive energy into the region. The fact that 
we see this behavior in the vacuum state as well suggests that it may represent the 
generic case. The correlations already seem to be "coded" in the vacuum state. 



A. Two Dimensions 



We begin with the calculation of the energy density-flux correlation function for a 
massless scalar field in the Minkowski vacuum state in two dimensions. The relevant 
components of the stress-tensor are 



T = ~T xt = - 



•ixVit J 



and 



T tt = T tt = l -[{ ( f )jt f + 



(45) 
(46) 

(47) 



Renormalized operators, such as : {(f), t) ', can be written as 

:(0,*) 2 :=(0,O 2 -((0,*) 2 )o 

where () denotes the vacuum expectation value (VEV). 

First use Wick's theorem to write the operator <f>, x <f>, t : (4>,t> ) : as a sum of normal- 
ordered products and VEVs: 



<P,x fat ■ (fat f ■ = ■ 4>,x fat (fat ) 2 : +2 : 4>, t 4>, t > ■ (4>, x <i>,t> > 

+ 2 : 4>, x (p, t , : (4>,t <f), v ) + 2(4>, x <f>, t , ) (<p, t <j>, v ) 



where we have also used the fact that <f), x (f>, t =: <fi, x <fi,t '■■ We can get (p, x cj), t : ((f), x > 
by letting t' — > x'. 

The flux-energy density correlation function in an arbitrary state is given by 
C = (: T xt (x) :: T*'*V) ■) = (T xt (x) : T*'*V) 
= ~2 (fa* fat ■ (fa* ) 2 + fax fat ■ {fax' ) 2 :> 



(48) 

2 
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(: <f), x <f), t {(j>, t , f :) + (: <f), x <f), t ((f), x > f :) 



+ 2(: <f),t (f), t > :) (4>, x <t>,t> ) + 2 (: <f>,t <f>,x> ) 
+ 2(: <f>, x <f), t > :) (0, t ) + 2(: 

ja:' )o 

+ 2{<f>, x <f), t > ) o (0, f ) + 2{(j), x (f), x > ) Q (<f>, t (j), x > ) t 



(49) 



Now consider the case where the given state is the vacuum. Then all of the VEVs of 
normal-ordered products vanish and we have 



(T x \x) : T 1 (x') :) 



■ <t>,t' ) (4>,t <P,t> ) + (<t>,x <P,x> ) O (0,t <fi,x' ) 



OJ ' 



with 



,x <j>,t> ) = d x d t iG = d x >d t G = (<p, t 4>, x > ) 



(50) 



(51) 



where G is the vacuum two-point function in two-dimensional Minkowski space. 
Therefore our correlation function becomes 



(T xt (x) : T 1 (x') :) = -(^ <j>, v ) ((0, t 0, t , ) + (0 j!e ) ) 
In two-dimensional Minkowski space we have that G = G(a) = 



-(4 



7T 



(52) 



In cr, 



where a = (l/2)(x" - x"*)^ = (1/2) [(a; - x') -(t - t') % and 9 x G(tr) = (x-x')G"((t) 
where G' = dG/da. Therefore we have that 

(favfa ) = d x d t ,G{a) = (x- x')(t - t')G"(a) 
{M,«) = d t ,[-(t-t')G\o)] = G\o)-(t-t'fG"(o) 



x0,,')o = d x ,[(x-x')G'(a)} = -G'(a)-(x-x') 2 G"(a) 



and so 



(<f>, t (j>, t , > + ((f), x (/), x , ) = -[(x- x'f + (t- t'f] G"{a) . 
The flux-energy density correlation function in the vacuum then becomes 

C = (T xt (x) : T*'*V) :> = (x - x')(t - t')[(x - x'f + (t - t'f] [G"(a)f 



(53) 
(54) 

(55) 



Let t — and x' = 0, i.e., we look at the correlation between the flux on the 
t — surface and the energy density on the x' — worldline. With these choices, our 
correlation function becomes 



C = -xt'[x 2 + r] [G"(a)Y 



(56) 



The sign of C in different regions of the x, t' plane are shown in Fig. [5j 

The correlation of the signs of the energy density, p, and flux, F = T xt , is shown 
in different regions of the x, t' plane in Fig. 
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t' 



C>0 



C < 



X 



C<0 



C>0 



FIG. 5: The sign of C in two dimensions is shown in different regions of the x, t' plane. 

Case A. For t' < and x > 0, then C > and for t' < and x < 0, then 
C < 0. Therefore in Fig. EJ^a), we see that if p > initially, then F{x < 0) < 
and F(x > 0) > subsequently, on average. In Fig. Mb), if p < initially, then 
F(x < 0) > and F(x > 0) < subsequently, on average. 

Case B. For t' > and x > 0, then C < and for t' > and x < 0, then C > 0. 
Therefore in Fig. [6](c), we see that if F(x < 0) > and F{x > 0) < initially, then 
p > subsequently, on average. In Fig. [6]^d), if F(x < 0) < and F(x > 0) > 
initially, then p < subsequently, on average. 

Thus in all cases, the sign of the energy density at x' = and the fluxes at t = are 
correlated as we would expect. 



B. Four Dimensions 

We now compute the vacuum energy density-flux correlation function in four- 
dimensional Minkowski space. The radial flux is given by 



F = T 



irt 



= W = 



-T, 



(57) 



and the energy density by 



P = t u = ~ {ci>, t <t>, t r). 



(58) 
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p>0 



(a) 



X 



t' 



X 



p <0 



(b) 



t' 

(i p>0 



(c) 



t' 



p <0 



x 



(d) 



FIG. 6: The correlation of the signs of the energy density and flux in two dimensions is 
shown in different regions of the x, if plane. 

For the correlation function we have 

C={: T rt (x) :: (x') :) = (T rt (x) : T*'*'(x') :) . (59) 
A calculation analogous to the one given for two dimensions yields 

C = --{: <f>,r<f>,t ■■■■ 0,t 0,* :) 

= -{<f>,r<f>,t' ) O (0>t0>f ) - (<P,t <!>,%' ) ( < l ) ,tV')) ■ (60) 

Here the vacuum two-point function is (4>(x)(f)(x')) = G(cr) = l/(87r 2 <r), where o = 
(1/2) [-At 2 + |Ax| 2 ], and At = t - t', Ax = x - x'. So we have that 

(4>, t 4>, t , ) = d t d t ,G{a) = a, tt , G'{a) + G"(a)a, t a, t , = G'(a) - (At) 2 G"(a) , (61) 
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and 



V0,t' ) = d r d t 'G(a) = a, rt > G'(a) + G"(a)a, r a, t > = Ata, r G"(a) 



Furthermore, 



\<t>,r<l>ri> )o 



-At a' 1 ' G"(a), 

a,i' r G'(a) + a,i> a, r G"(a) 



so 



(0, r fa ) o (0, t 0'*')o = " A * ^"(a) (a' 1 ' a nV ) G'(a) + (a' 4 ' cv ) a rr G"(a) 



(62) 



(63) 



(64) 



We need to calculate the derivatives of a which appear in the above expression. 
The scalar quantity a' 1 cv is most easily computed in Cartesian coordinates, where 
a„i = —Ax* and hence a' 1 a„i = |Ax| 2 . We consider the case when x' = 0. We then 



find 



<T, r 



CP <T,i' 



r 
r 



As a result our energy density-flux vacuum correlation function reduces to 



C = r At 



(At) 2 + r 2 [G"(a)Y 



(65) 



(66) 



where a = (1/2) [r 2 — (At) 2 ] for r' = 0. Note that the sign of C depends only on the 
sign of At. Recall that At — t — t'. If t > t', then C > 0. Consequently, p(t') > 
implies F = T rt (t) > on average, as illustrated in Fig. [7J and p(t') < implies 
F = T rt (t) < on average. This is consistent with our two-dimensional results. If 
p > at the origin at one time, the flux at a later time is likely to be outgoing. 
Similarly, if p < at one time, the later flux is more likely to be ingoing. We can 
also reverse the time order. An outgoing flux at t is correlated with p < at t' > t, 
and an ingoing flux at t is correlated with p > at if > t. 

Lastly we discuss the falloff of the correlations in both two and four dimensions. 
In two dimensions, choosing t — 0, we have At — t — if — —if and so 



Since in two dimensions G = —(Att 
let x — > oo, with t' fixed, then 



C=-xt' (x 2 + t' 2 ) [G"(a)] 2 . 

1 In a, G" = l/(47rcr 2 ) oc -l/(a; 2 

t' 



c ~ — 



X'' 



Alternatively, if we let t' — > oo, with x fixed, then 

x 



C 



5 ' 



(67) 
t' 2 ) 2 . If we 

(68) 
(69) 
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FIG. 7: Here the relation between the energy density at time t' and r' = and the flux F at 
time t and radius r is sketched. A positive energy density at an earlier time t' is correlated 
with an outgoing flux at a later time t, and vice versa. 

Thus the falloff is symmetrical for both spacelike and timelike separations. 

In four dimensions, with r' — and since G = l/(87r 2 cr), G" = l/(47r 2 cr 3 ) oc 

l/[r 2 — (At) 2 ] , we have that for r — > oo, with At fixed, then 

C-^- (70) 



r 9 



For At — > oo, with r fixed, 



C q. (71) 

(At) 9 V ; 

and again the falloff is symmetric between space and time. Note also that C has the 
same form in two and four dimensions, apart from the difference in the expression 
forG(tr). 



IV. CONCLUSIONS 

In this paper we have examined in detail the behavior and correlation of energy 
density and flux in a "Heifer state". This state, first suggested by Adam Heifer, is 
one in which the energy density can be made arbitrarily large and negative over an 
arbitrarily large region of space, even in Minkowski spacetime. However, this can be 
achieved only at one instant in time. For the quantum inequalities to hold, there 
must be large time-varying fluxes associated with these states. We show that this 
is in fact the case, and analyze the correlated behavior of the expectation values of 
the flux and energy density, in order to gain a deeper insight into these remarkable 
quantum states. 
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In the second part of the paper, we examine the energy density-flux correlation 
function in the Minkowski vacuum state. Although we now work with correlation 
functions, as opposed to expectation values, we find qualitatively similar behavior for 
energy density and flux fluctuations. An initial negative energy vacuum fluctuation 
in some region of space is correlated with a subsequent flux fluctuation of positive 
energy into the region. 

The quantum inequalities place a (negative) lower bound on the distribution func- 
tion for energy density vacuum fluctuations |26|. Perhaps the mechanism which 
ensures that the quantum inequalities hold in the Heifer state, as well as in other 
quantum states associated with negative energy, is, at least in some sense, already 
"encoded" in the fluctuations of the vacuum. 
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